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We study numerically the properties of local low-energy excitations in the two-dimensional Ising 
spin glass. Given the ground state, we determine the lowest-lying connected cluster of flipped spins 
containing one given spin, either with a fixed volume, or with a volume constrained to lie in a certain 
range. Our aim is to understand corrections to the scaling predicted by the droplet picture of spin 
glasses and to resolve contradictory results reported in the literature for the stiffness exponent. We 
find no clear trace of corrections to scaling, and the obtained stiffness exponent is in relatively good 
agreement with standard domain wall calculations. 
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I. INTRODUCTION 

According to one of the principal scenarios for the spin 
glass state, the "droplet picture"i242i4, the minimum en- 
ergy excitation of hnear dimension L containing a given 
spin, a droplet, has an energy proportional to . The 
"stiffness exponent" 6 is positive if one assumes that the 
transition temperature Tc is finite. It follows that, in the 
thermodynamic limit, excitations that flip a finite frac- 
tion of the spins cost an infinite energy. However, the re- 
sults of several numerical calculations^iSiSiSiifliiiiiSi, on 
small system sizes, imply that the amount of energy 
needed to generate system-size droplet excitations is in- 
dependent of size, so that 9 for droplet-like excitations 
is zero. This is consistent with the alternative replica 
symmetry breaking (RSB)i2iiiiiSii& picture. Note that 
numerics simultaneously finds that the stiffness exponent 
for domain-wall excitations is positive and therefore ap- 
parently different from the exponent for droplets, an un- 
expected feature. 

To explain the discrepancy between those two dif- 
ferent estimates of the stiffness exponent, it has been 
proposediii^ that the droplet theory is correct on large 
scales, i.e. the stiffness exponent 6 is the same for droplets 
and domain walls, and that the apparently contradictory 
numerical data can be e xplaine d by corrections to scal- 
ing. More precisely Refs. llTllSl propose that the droplet 
energy AE has the form 

AE = AL^ + BL-'^, (1) 

where 9 > {if Tc > which is the case for d > 3), 
w is a correction to scaling exponent, and A and B are 
positive constants. If > 0, AE has a minimum at 
some value of L and if this occurs in the range of sizes 
where numerical data is taken, the numerical values could 
be roughly independent of size in the range studied. In 
general, only a modest range of sizes, L < 10, can be 
studied in dimension d equal to 3. 



To quantitatively check this proposal, simulations have 
recently been carried out^® on a one-dimensional model 
with long range interactions which fall off as a power of 
the distance, in a region of parameters for which Tc > 
0. A much larger range of sizes can be studied than in 
three dimensions but even up to L = 512 the energy of 
system-size excitations was found to be independent of 
size. Hence, for Eq. |^ to be correct with 9 > the 
values oi 9, u, A and B must conspire to give an almost 
constant AE over a very wide range of sizes. 

Another model where large sizes can be studied is the 
two-dimensional Ising spin glass at T = 0, where effi- 
cient algorithms^fiiSiiSS can be applied to determine exact 
ground states. Now 9 < in d = 2, with T = "domain- 
wall" calculations giving^^i^3,24,25,26,27 q ^ _q 23. How- 
ever, finite-T simulationsSSiSS and T — studies in a 
magnetic field'^*', both of which excite droplets rather 
than domain walls, find 9 ~ —0.47. A good review of 
this situation is given in Ref. 31. Since 9 < in d = 2, 
both the and L~'^ contributions to the dropl et ene rgy 
in Eq. Q decrease with increasing L, but Refs. ar- 
gue that they could combine to give an effective exponent 
of about —0.47 for small L crossing over to the asymp- 
totic value of about —0.28 for larger sizes. Hartmann and 
MoorSiiS find evidence for this crossover in their numer- 
ical studies in which a certain prescription was used to 
generate droplets, though earlier Hartmann and one of 
ua^, who generated droplets in a different way, did not 
find evidence for an effective exponent close to —0.47. 
This difference may indicate that the amplitude of the 
correction term B in Eq. depends significantly on the 
precise way in which the droplets are generated. In yet 
another way of generating droplets, Picco et ali^^^ find 
9 close to the —0.47 value but with no crossover to the 
domain wall value of around —0.28 for the larger sizes. 

Recently, in an interesting paper, Lamarcqet al^ have 
calculated the energy of droplet (cluster)^ excitations in 
three dimensions, in a manner very similar to the spirit 
of the droplet mode l-'- Interestingly they find that the 
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energy actually decreases, although very slightly, with 
increasing size. They also find that the volume of the 
droplets has a non trivial fractal dimension less than the 
space dimension. In Ref. fl^i this slow decrease of the 
droplet energy was taken as a possible evidence of the rel- 
evance of the correction to scaling term in Eq. Note 
however that no sign of a crossover towards the suppos- 
edly correct value of the stiffness exponent was reported, 
even when larger clusters could later be included in the 
analysis^. 

Here we perform a calculation similar to that of 
Lamarcq et al. but in two dimensions. The aim is to 
see if the energies of the droplet energies calculated in 
this way fit Eq. |^ and give an effective exponent close 
to —0.47 at small sizes crossing over to about —0.28 at 
larger sizes as has been proposed by Mooreii and Hart- 
mann and Moored. Since the way our clusters are gen- 
erated is directly inspired by the original definition of the 
droplets, one could expect the prediction to be well 
suited in our case. Instead, we do not find any crossover 
in AE by applying the definition of Lamarcq et al. in 2d 
for clusters as large as 64 spins, neither do we find a very 
negative exponent —0.47. We find however an exponent 
similar to —0.47 at very small sizes if a slight modifica- 
tion of the definition of the clusters is made, namely if the 
size of the clusters is let free to evolve between n and 2?!, 
as suggested by the original "scale invariant" definition 
of the droplets. 

II. THE MODEL AND DETAILS OF THE 
NUMERICS 

We take the standard Edwards-Anderson spin glass 
model 

H ^ -^Ji-jSiSj, (2) 

where the Si — ±1 are Ising spin variables at the sites 
of a simple cubic lattice, and the Jij are nearest neigh- 
bor interactions with a Gaussian distribution with zero 
mean and standard deviation unity. Periodic boundary 
conditions are applied on lattices with N — LF' spins. For 
most of our work we take L = 64 but wc also did some 
calculations with smaller sizes down to L = 16. We use 
Ng = 1000 realizations of the disorder. 

We are interested in the properties of low-energy, 
droplet-like excitations above the ground state of the 
Hamiltonian in Eq. |2Jl. In order to determine these, 
we first find the ground state for each realization of the 
disorder. This is done using either the Koln spin glass 
server—, or, for small sizes, by parallel tempering'^ME. 

Following Ref. |3^ we then generate a droplet by first 
choosing randomly a "central" spin in the system and 
reversing it. We then construct a connected cluster of 
fixed size n around this central spin by fiipping all the 
spins in this cluster. We then find the new ground state 
with the following three constraints: 
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FIG. 1; Randomly chosen clusters of size n — 64 obtained by 
the procedure described in Sec. ITTl 

1. The central spin is always fiipped with respect to 
the ground state, 

2. The number of spins in the cluster is always n, 

3. The cluster is always connected. 

The new ground state is found by a combination of par- 
allel tempering and a Kawasaki-type dynamics for the 
spins in the cluster, in order to conserve its size n con- 
stant. Our algorithm does not satisfy detailed balance, 
but this is irrelevant since we only want to find the new, 
constrained ground state. The lowest temperature in 
the parallel tempering is taken to be sufficiently low, 
T = 0.02, that the system is always in a local minimum 
with no random noise kicking the system to higher en- 
ergy states. Our results for ground states in the presence 
of a droplet are obtained from spin configurations at this 
lowest temperature. We study cluster sizes up to n = 64. 

In this procedure, which follows Ref. IssL the size of 
the cluster is strictly fixed. It is therefore slightly dif- 
ferent from the droplets of Fisher and Hus o^i^i^ which 
are defined in a "scale-invariant" way, i.e. the size is not 
strictly fixed but allowed to vary over a certain scale. 
We have therefore also computed the energies of scale- 
invariant excitations in which the size of the cluster is 
allowed to vary by a factor of 2, more precisely to lie be- 
tween n/2 and n — 1 for n = 4, 8, • ■ •, 64. Examples of 
some randomly chosen clusters of size n — are shown 
in Fig.Hl 

Our criterion for ensuring that we have found the true 
ground state in the presence of an overturned cluster is 
as follows. We first make a few (typically 10) extremely 
long runs with different bond configurations to estimate 
the typical time scale to reach the new ground state. 
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^typ {n) . We record the energy versus time of the lowest 
temperature for the two copies, for simulations of at least 
lOttyp sweeps. From these graphs, we can evaluate ^typ 
as the typical time it takes those 10 realizations of the 
disorder to give a result that can be confidently taken 
as the ground state. Then, for each of the Ng samples, 
we first run two copies of the system with independently 
drawn initial clusters for ttyp{n) Monte Carlo steps. We 
then require that the two copies have the same energy 
for a further continuous period of time equal to ttyp{n). 
In this way, we spend more time on the "hard" samples, 
which need several times ityp (n-) to converge, than on the 
easy samples for which we do not need to run much more 
than ttyp{n) sweeps. Obviously, one can not be absolutely 
sure that the ground state has really been found since it 
might happen that both copies stay for a time larger than 
^typ("-) in the same state which is not the ground state. 
However, this seems to be fairly unlikely. To be on the 
safe side we used ttyp{n = 64), for all sizes, even though, 
for n < 64, one has ttyp{n) <^ ttyp{n — 64). 

For the scale-invariant clusters we start one copy with 
the minimum-size cluster, n/2, and the other with the 
maximum-size cluster, n — 1, to minimize further the 
probability that the two copies inadvertently spend a 
long time in the same state which is not the ground state. 

We now present our numerical results. 

III. RESULTS 
A. Clusters of Fixed Size 

In this section we present results of the simulations 
where we fix the droplet size precisely, as was done by 
Lamarcq et al.-?- on the 3d model. 

1. Droplet Energies 

Fig. 121 shows the droplet energy, {AE), where the av- 
erage is over Ng = 1000 samples, as a function of n, the 
number of spins in the droplet, on a log-log scale. For 
each sample, AE is the energy difference between the 
constrained and unconstrained systems, by definition a 
positive quantity. If we define the fractal dimension of 
the volume of the droplets to be df, i.e. 

iir^R^f, (3) 

where R is the mean "radius" of the droplet (to be defined 
below) then the slope of the data is fl/d/, since by defini- 
tion {AE) ^ . The expectation of the droplet theory 
is that df is equal to the space dimension, df = d = 2. 
The slope of -0.115 for i = 64 would lead to 6 = -0.23, 
not quite as negative as the established value of —0.28 
from domain wall calculation, but perhaps the difference 
is not significant given the rather small range of droplet 
sizes than we can study. Actually, fits to data presented 
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FIG. 2: Mean droplet energy against the number of spins in 
the droplet for lattice sizes L — 16 and 64 for droplets of fixed 
size. Data for L = 32 is very similar to that for L — 64 and 
is not shown for clarity. The slope for L = 64, which is equal 
to 9/df, is —0.115. No crossover in the slope can be detected 
in the L = 64 data. 



below indicate that d/ is somewhat less than two, which 
would accentuate the difference between our value of 9 
and —0.28. However, the apparent difference between df 
and 2 may itself be due to corrections to scaling. 

We note also in Fig. |21 that the data for L = 16 shows 
an interesting feature. For small cluster sizes, the rela- 
tion between (AE) and n follows that of the L = 64 
and the two sets of data start to depart from each other 
at large cluster sizes, so that the apparent stiffness expo- 
nent becomes more negative at large sizes. This crossover 
is however exactly opposite to the one expected from 
Eq. . A visual inspection reveals that even the clusters 
with n ~ 64: have a radius smaller than the system size 
L = 16, so that the effect is not trivially due to the fact 
that boundaries of the clusters interact with each other. 
Rather, this crossover is most probably due a finite size 
effect on the initial ground state energy, a smaller size 
making it too large, resulting then in too small an en- 
ergy difference AE. 

An important conclusion drawn from Fig. [21 is that no 
crossover from a large to a smaller value of the stiffness 
exponent is visible, although our data cover a reasonably 
large size window. We are led to the conclusion that 
the definition of the droplets used here is relatively free 
of the corrections to scaling expected from Eq. and 
moreover that it leads, for d = 2, to a determination of 
the stiffness exponent in fair agreement with domain wall 
calculations. 
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FIG. 3: Distribution of the energy of the droplets of fixed size 
n, for different n and system size L — 64. 
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FIG. 4: The data in Fig.^scaled by the mean energy (AE), 
according to Eq. 



2. Distribution of Droplet Energies 

Fig. O shows the distribution of the energy, P{AE), 
of fixed size droplets for different sizes with L = 64. In 
Fig. ^ we show that these distributions can be satisfac- 
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FIG. 5: The number of spins in the droplet, n, as a function 
of various definitions of the "radius", as described in the text. 
This is for droplets of a fixed size. 



torily scaled using the form 



(4) 



using the mean energy (AE) which is plotted in Fig. |2 
The data indicates that the scaling function P{x) varies 
linearly for small x. This result was also obtained in 
three-dimensions^^, but differs from the usual assump- 
tion in droplet theory that the distribution has a finite 
weight at the origin. 



3. Fractal Dimension of the Droplets 

To determine the fractal dimension of the volume of 
the droplets through Eq. (O , we need to give a definition 
of their "radius" . We have done this in different ways and 
checked that they give consistent results for df. We mea- 
sured Rx and Ry, the maximum extent of the droplets in 
the X and y direction as well as max(i?2;, Ry). In addition 
we determined the radius of gyration, Rg, of the droplet 
defined by: 



Rn 



n ^ 

itE cluster 



^ cm I J 



(5) 



where rem is the position of the center of mass of the 
cluster. 

Figure 13 shows the result of this analysis for L = 64, 
together with power law fits to the data. The agreement 
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FIG. 6; The mean energy of the scale-invariant droplets and 
system sizes L = 16 and L — 64. For the size indicated as n, 
the droplets are allowed to range in size between n — 1 and 
n/2. 

between Rg and T[iax{Rx,Ry) is excellent. If we exclude 
the L = 4 point, the fits give df closer to 2. For example 
the fit to Rg gives df ~ 1.88. Combining these values of 
df with 0/df obtained from Fig. [3 we get: 

{—0.23 (assuming df — 2), 
-0.21 (fit to Rg), (6) 
-0.22 (fit to Rg, excluding n = 4). 

We emphasize again that no trace of a more negative ex- 
ponent of the order of —0.47 can be detected in our data, 
thus in disagreement with the predictions by Mooreii 
that this method should reveal the corrections to scaling 
included in Eq. 

B. Clusters with a Range of Size 

In this section we allow the size of the droplets to vary 
over a factor of two, more precisely between n/2 and 
n — 1 . This is more in the spirit of the droplet picture of 
Fisher and Huse^^, where droplets are defined, in 3d, 
as objects with boundaries lying between a cube of linear 
size £ and another one of size 2£. We now perform the 
same analysis as above for these new clusters. 

1. Droplet Energies 

The mean droplet energy as a function of size n is 
shown in Fig. El For the size specified as n, the droplets 
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FIG. 7: Distribution of the energy of droplets of variable size, 
for different ranges of size for L = 64. For a size indicated as 
n, the droplets are allowed to range in size between n — 1 and 
n/2. 



were allowed to vary in size between n — 1 and n/2. 
For very small sizes, where the data for L = 16 and 
64 agree, the slope is about —0.22, considerably more 
negative than the value of —0.115 found for droplets of 
fixed size in Sec. IIII Al For larger sizes, where the data 
for i = 16 differs from that for L = 64, the slope of the 
presumably more reliable i = 64 data becomes less neg- 
ative and is compatible with —0.115, though there is not 
enough range of size to allow for a precise fitted value. If 
we insert df = 2, then, according to Eq. this would 
give a crossover from 9 ~ —0.44 to 9 ^ —0.23, not very 
different from the crossover proposed in Refs. fT^ITM 

However, it should be pointed out that no such 
crossover was found in Sec. IIII Al for fixed size droplets, 
and, furthermore, the values of the droplet radius R, dis- 
cussed in Sec. IIII B 31 are extremely small, of order one 
or two lattice spacings, in the region where the data of 
Fig. suggests 6 ~ —0.44, so that the meaning of a 
power law correction of the type for such sizes is 
not clear. Note finally that the apparent fractal dimen- 
sion for these very small clusters is quite far from 2 (see 
below), so that, if we accept this value for df the initial 
apparent stiffness exponent is around —0.35, rather far 
from the value —0.47 expected by Mooreiii^. 

Finally, we have to admit that we do not understand 
why the crossover in Fig. is absent when the size is 
exactly fixed, see Fig. (21 
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FIG. 8; The data in Fig. Qscaled by the mean energy (AS), 
according to Eq. (0J. 



2. Distribution of Droplet Energies 



FIG. 9: The number of spins in the droplet, n, as a function 
of various definitions of the "radius", as described in the text. 
This is for droplets whose size is allowed to vary between n—1 
and n/2. Data for n — 4 has been excluded from the power 
law fit. 



The distribution of droplet energies is shown in Fig. d 
and the scaled data is plotted in Fig. |S1 again making use 
of Eq. 0). As for droplets of fixed size, the distribution 
is essentially linear for small energy, though in this case 
there is a very small, but definitely non-zero, intercept. 



3. Fractal Dimension of the Droplets 

Figure |51 shows the various definitions of the radius of 
the droplets, discussed in Sec. 1111 ^TSl for different values 
of n. The data does not fit well a straight line if the 
smallest size is included (n = 4, corresponding to droplets 
with 2 and 3 spins) so this has been omitted in the fits. 
The values of c?/ for Rg and max(_R2;, Ry) agree well, as 
for droplets of fixed size, though the value of df (~ 1.7) is 
a little lower than in that case. However, the trend is for 
the slope to increase with increasing n, this is particularly 
noticeable for the Rg data, and so it is possible that df 
is actually equal to 2, as we concluded also in the case of 
the cluster of fixed size. 



IV. CONCLUSIONS 

The main motivation for our work is to test Moore's 
proposalii that contradictory results of the literature, 
both in 2 and 3 dimensions, for the determination of the 
stiffness exponent 9 in spin glasses were mainly due to 



"corrections to scaling" and captured by Eq. (Q. We 
have performed simulations similar to those of Lamarcq 
et al.?'^-'^^ but in the 2d case, where this theoretical con- 
jecture could be checked. From our results, we can draw 
the following conclusions. 

1. Following exactly the procedure of Ref. Issl we do 
not find any crossover in the behavior of the energy 
of the cluster as a function of their size, up to n = 
64, contrary to Moore's prediction; 

2. The value of the extracted stiffness exponent, 9 ~ 
—0.23, is in fair agreement with domain wall calcu- 
lations, 6 — —0.28. This result is clearly different 
from the 3d case, where the droplet and domain 
wall exponents are found to be very different. 

3. For scale invariant clusters, a crossover not incon- 
sistent with ^ is obtained, but only for very small 
cluster sizes. 

As a conclusion, the problems we wanted to tackle in 
this paper unfortunately remain unsolved. We still do 
not understand the 3d results by Lamarcq et al.'^^ and we 
do not understand the origin of the two different stiffness 
exponents reported in 2d. Apart from very small sizes for 
one of the two methods of generating clusters, we did not 
find evidence for the crossover implied by Eq. (Q. The 
nature of corrections to scaling in spin glasses therefore 
remains poorly understood. 
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Suppose we take the point of view that the differing 
claims obtained in Refs . Il^32ll35i34 and the present work 
for the existence, or otherwise, of the proposed crossover 
in two-dimensions, are evidence for a wide range of values 
of the amplitude of the correction B in Eq. It is then 
surprising that, to our knowledge, all estimates of 9 from 
Monte Carlo simulations at T <C Tc give 9 (for droplets) 
close to (and consistent with) zero, and different from 
9 for domain walls which is positive. This is true for 
Ising spin glasses in three and four dimensions^", vector 
spin glass model a^i^" , and a one-dimensional model with 
long range interactionsiS. In addition, various T = 
calculations^iSii also find 9 (droplets) ~ 0. If the ampli- 
tude of the correction term is highly non-universal, as 
perhaps implied by the 2d results, it is surprising to us 
that the parameters in Eq. systematically conspire to 
give a droplet energy for, say, Monte Carlo simulations 
which is independent of size over the range studied. Note 
too, that for the Id simulations^^, the range of sizes is 
very large, up to i = 512. 

Of course, the usual criticism that the truth lies be- 
yond the reach of numerical simulations may be true. In 
this case, it is worth discussing the relevance of physics 
which can only be seen on such astronomically long time 
scales. From a mathematical point of view, understand- 
ing this physics is of fundamental importance. However, 
if the asymptotic behavior would only be seen on time 
scales which are inaccessible to experiments as well as to 
simulations, then the physical relevance is rather limited. 
In fact, experimental systems below the spin glass transi- 
tion temperature Tc are not truly in equilibrium but are 
only equilibrated up to some coherence length scale, £{t), 
which increases slowly, presumably logarithmically, with 



the time of the experiment t and also depends on T. One 
of us and Bouchaud^ have estimated that £{T) < 10 (in 
units of the spacing between spins) at reasonable exper- 
imental time scales at T ~ 0.5Tc. Refs. 0!43! estimate a 
somewhat larger value (rather more than 20) at this tem- 
perature and quite a lot larger^^ closer to Tc- Hence a 
reasonable estimate of £{t) well below Tc in experiments 
is around 10-20 times the spacing between the spins, not 
very different from the sizes that can be simulated nu- 
merically, and much less than the sizes which appear to 
be necessary to see droplet behavior. 

It is therefore possible that even if the droplet theory 
is asymptotically correct, the region in which its predic- 
tions can be observed quantitatively is not accessible ex- 
perimentally. Close to Tc the length scales which can be 
equilibrated in experiments may be much larges^ii^, but 
in this region critical fluctuations give significant correc- 
tions to droplet behavior so even larger length scales are 
needed^*^ to see droplet behavior than at lower temper- 
atures. 
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